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Abstract. This article deals with the study of Bianchi type-I universe in the context of f(R, T) gravity. 
Einstein’s field equations in f(R,T) gravity has been solved in presence of cosmological constant A and 
quadratic equation of state (EoS) p = ap 2 — p, where a ^ 0 is a constant. Here we have discussed two 
classes of f(R : T) gravity i.e. f(R : T) = R + 2 f(T) and f(R,T) = fi(R) + h(T). A set of models has 
been taken into consideration based on the plausible relation. Also we have studied the some physical and 
kinematical properties of the models. 
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1 Introduction 


It is known that [lj|5] in present scenario, our universe is accelerating. However, final satisfactory explanation about 
physical mechanism and driving force of accelerated expansion of the universe is yet to achieve as human mind has 
not achieved perfection. From the modern cosmology, it is known that a point of universe is filled with dark energy. 
It has been addressed by various slow rolling scalar fields. It is supposed that the dark energy is responsible to 
produce sufficient acceleration in late time of evolution of the universe. Thus it is much more essential to study the 
fundamental nature of the dark energy and several approaches have been made to understand it. The cosmological 
constant is assumed to be the simplest candidate of dark energy. It is the classical correction made to the Einsteins 
field equation by adding cosmological constant to the field equations. The introduction of cosmological constant to 
Einsteins field equation is the most efficient way of generating accelerated expansion but it faces serious problems 
like fine-tuning and cosmic coincidence problem in cosmology [6,7. Quintessence [8], phantom [9], k-essence 10 


tachyons [TTJ , Chaplygin gas [12] are the other representative of dark energy. However, there is no direct detection of 
such exotic fluids. Researchers are taking interest to explore dark energy due to lack ofstrong evidence of existence of 
dark energy. Several authors (Singh et al. [l3], Singh et al. 
model with cosmological constant in different contexts. 

Dark energy can be explored in several ways and modifying the geometric part of the Einstein-Hilbert action 16 


14 , Jamil and Debnath 1151) have discussed cosmological 


is treated as the most efficient possible way. Based on its modifications, several alternative theories of gravity came 
into existence. Some of the modified theories of gravity are f(R ), /(T), /(G) and f(R,T) gravity. These models are 
proposed to explore the dark energy and other cosmological problems. Sharif and Azeem 17 discussed the Cosmological 
evolution for Dark Energy Models in f(T) Gravity. The f(R) modified theory produces both cosmic inflation and mimic 


behavior of DE including present cosmic acceleration 18 -20 . Amendola et al. 21 have discussed the cosmologically 


viable conditions in f(R) theory which describe the dark energy models. Recently author like Sahoo and Sivakumar 22 
Ahmed and Pradhan 231, Pradhan et al. 24 have investigated the cosmological models with cosmological constant 


in /(i7, T) gravity for different Bianchi type space-time. 
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Quadratic equation of state is needed to explore in cosmological models due to its importance in brane world model 
and the study of dark energy and general relativistic dynamics for different models. The general form of the quadratic 
equation of state is given by 

p = p 0 + ap + /3p 2 , (1) 


where po,a,/3 are parameters. This equation 0 is nothing but first term of Taylor expansion of any equation of state 
of the form p = p(p) about p = 0. 

Ananda and Bruni discussed the cosmological models by considering different form of non-linear quadratic equation 
of state. Ananda and Bruni || 25] have investigated the general relativistic dynamics of RW models with a non-linear 
quadratic equation of state and analysed that the behaviour of the anisotropy at the singularity found in the brane 
scenario can be recreated in the general relativistic context by considering a equation of state of the form 0 . Again 


they 26 have discussed the anisotropic homogeneous and inhomogeneous cosmological models in general relativity 


with the equation of state of the form 


p = ap -f-, 

Pc 


( 2 ) 


and tried to isotropize the universe at early times when the initial singularity is approached. In our present study, we 
have consider the quadratic equation of state of the form 


p = ap - p, 


(3) 


where a ^ o is a constant quantity and such type of consideration does not effect the quadratic nature of equation of 
state. 

Nojiri and Odintsov 


27 studied the effect of modification of general equation of state of dark energy ideal fluid by 


the insertion of inhomogeneous, Hubble parameter dependent term in the late-time universe. The quadratic equation 
of state may describe the dark energy or unified dark energy (27,128]. Rahman et al. 29 investigated the construction 


of an electromagnetic mass model using quadratic equation of state in the context of general theory of relativity. Feroze 
and Siddiqui 30 studied the general situation of a compact relativistic body by taking a quadratic equation of state 
for the matter distribution. Maharaj and Takisa 31 have investigated the regular models with quadratic equation of 


state. They have considered static and spherically symmetric space-time with a charged matter distribution and found 
new exact solutions to the Einstein-Maxwell system of equations which are physically reasonable. 

A Cosmological model based on a quadratic equation of State unifying vacuum energy, radiation, and dark energy 
has been discussed by Chavanis 32 and also a cosmological model describing the early inflation, the intermediate 


decelerating expansion, and the late accelerating expansion by a quadratic equation of state have been investigated 
by the same author 33 . Strange Quark Star Model with Quadratic Equation of State has been investigated by 
Malaver 34 and they have obtained a class of models with quadratic equation of state for the radial pressure that 
correspond to anisotropic compact sphere, where the gravitational potential Z depends on an adjustable parameter 
n. Recently Reddy et al. 35 has studied the Bianchi type-I cosmological model with quadratic equation of state in 


the context of general theory of relativity. 

Motivated by the aforesaid research, we have investigated the Bianchi type-I cosmological model in /(R, T) gravity 
with quadratic equation of state and cosmological constant . Here we have discussed two classes of /(R, T) gravity. 


2 Gravitational field equations of f(R,T) modified gravity theory 

Let us consider the action for the modified gravity as 

s =/(w + l ")^ a ' < 4 » 

where /(R, T) is the arbitrary function of R and T. R is the Ricci scalar and T is the tress of the stress-energy tensor 
of the matter T^-. L m is the matter lagrangian density. For the choice of /(R, T ) we will get the action for the different 
theories . If /(R, T) = f(R) and /(R, T) = R then Q represents the action for f(R) gravity and general relativity 
respectively. The stress energy tensor of matter is defined as 

rp _ _ 2 5{y/=gL m ) 

v=g 6g ij ’ ( j 

and its stress by T = g^Tij. If we consider that the matter lagrangian density L m of matter is depends only on gij 
and not on its derivatives leads us 

dL m 
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By varying the action 0 w.r.t the metric tensor component gp we have 

f R (R, T)Rij - T) 9ij + ( gij D - ViVj)f R (R, T) = 8 tt - fr(R, T)Ty - MR, T)G ij , 


( 7 ) 


where 


@ij — —2Tij 


" 9 ijL n 


2 g lk 


O^Lrn 

dg^dg lk 


(8) 


Here /t(R,T) = Jr{R, T ) = , □ = V z Vj is the De Alembert’s operator and T; t j is the standard matter 

energy momentum tensor derived from the lagrangian L m . By contracting equation (|7|), we obtained the relation 
between R and T as 


fn(R, T)R + 3D f R (R, T ) - 2 f(R, T ) = SnT - f T (R , T)T - / T (ii, T)0, (9) 

where & = g l ^&ij. From 0 and the gravitational field equations can be written as 

f R (R, T)(Rij - l -R 9ij ) + ±f(R, T) gi j = (8tt - f T (R, T))(T id - l -T 9ij ) - f T (R , T)(©y - l©<7y) + O (10) 

The perfect fluid form of the stress energy tensor of the matter lagrangian is given by 


T ij = (p + p)^j -P9ij, (11) 

where iP = (1,0, 0,0) is the four-velocity vector and satisfies the relation idiq = 1 and u l X jUi = 0. p and p are the 
energy density and pressure of the fluid respectively. From equation ^ we have 

Oij = - 2 - pgij ( 12 ) 


It is to note that the functional f(R,T) is depends on the physical nature of the matter field through the tensor 0^. 
Thus for each choice of f(R,T) leads us to different cosmological models. Herko et al 36 presented three class of 
f(R,T) as follows 

R + 2/(T) 

f(R,T)=lf 1 (R) + f 2 (T) (13) 

Ji(R) + f2(R)fs(T) 

In this present work we have discussed two class of f(R,T) i.e. f(R,T ) = R + 2 f(T) and f(R,T) = fi(R) + / 2 (F). 
For the choice of f(R, T) = R + 2 f(T) and with the help of © and ( |l2| ), equation 0 takes the form 


Gij = (8tt + 2 f'(T))Tij + (2 pf(T) + f(T))g k 


(14) 


Which is the gravitational field equation in f(R,T ) modify gravity for the class f(R,T) = R + 2/(T). For the choice 
of f(R,T) = fi(R) + / 2 (F) and with the help of ( [lT] ) and (12), equation 0 takes the form 


f[(R)Ri 


pi(R)dij + (dij a - VjVj)/((i?) 


(Stt + f 2 (T))Tij + (/'(7> + -h{T))g, 


(15) 


Which is regarded as the gravitational field equation in f(R,T) modify gravity for the class f(R,T) = fi(R) + fiiT). 


3 Field equations and cosmological model for f(R, T) = R + 2 f(T) 

In f(R, T) theory, the gravitational field equations © in presence of Cosmological constant A are given as 

= [Sir + 2 f'(T)\ + [2 pf'(T) + f(T) + A] g ij: (16) 


where prime denotes differentiation with respect to the argument. For the choice of f(T) = AT ,eqn.(16) takes the 
form 

Gij — + 2A)] + [A p — pX + A] gp (17) 

Let us consider the Bianchi type-I space-time in the form 


ds z = dt z - X\dx 2 - X z dy z - X z dz z , 


2 7,2 


-2 r2 


(18) 
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where Xi, X 2 and X 3 are function of t only. The field equationsfor the line element (18) takes the form 


XiX 2 , XaX 3 , X2X3 , c , x „ 

( 87 r + 3X)p + pA A 

A1A2 A1A3 A2A3 

( 19 ) 

v + Y + Y Y — A p A 

A 2 A3 A2A3 

(20) 

Xi v 3 Vx 3 , 8 , , 

v + Y + Y Y — ( 8?r + 3 ^)P V ^ 

A1 A 3 AiA 3 

(21) 

x + x + x x — ( 8?r +V ^ 

A1 A 2 A1A2 

(22) 


4 Solution Procedure 


Now our problem is to solve the Einstein’s modified field equations (19)-(22). Here the system having four equations 
and six unknowns (Xi, X 2 , X 3 , p , p and A). To obtain the complete solution, we need two more physically plausible 
relations. The considered two physically plausible relations are 


1. Quadratic equation of state 

2 . Expansion Law 

— Power law 

V = V 0 t 3n , (23) 

— Exponential law 

V = V o e 0t , (24) 

where n and /3 are poitive constant quantity. According to the choice of expansion Law,we have obtained two different 
models of the Bianchi type-I universe. 


4.1 Power law Model 


With the help of equations (20)-(22),we have obtained the metric potentials as 

[/ 


Xi(t) = XioV 3 exp 


Xm 

V 


, * = 1,2,3, 


(25) 


where Wo and are constant of integration (i = 1,2,3) and which satisfies the relation Yl^=i Wo = 1 and Y^i=i = 
0. From the equations (19)-(20) and along with © we have got 

o 1 


a(87r + 2 A) 


X 2 X 3 X 2 X 3 X-, X 2 X 1 X 3 X 2 X 3 
X 2 + X 3 + X 2 X 3 X a X 2 XiXs x 2 x 2j 


Using (23) in (25), we have the metric potential as 

Xi(t) = X i0 V^exp 


—Xq it 


—3n+l 


(3 n - 1)V 0 


, i = 1,2,3 


(26) 


(27) 


The directional Hubble parameters are obtained as Hi = j + , i = 1, 2, 3. The Hubble parameter (iJ), deceleration 

parameter (q) , expansion scalar (0) and Shear scalar (a 2 ) are as follows 


^ _ 1 ^ 9 *ns + *na + *02*03 


H = —, q = —1 -\—, 0 = 3 — , cr = 
t n t 


V 2 t^ 


Using the observational value for q = —0.33 ± 0.17 37 , we have restricted n as n G (1.19,2) in case of Power law 
model. Here we noticed that H,0 and a 2 die out for larger values of t. With the help of ( |27| ) from (26), the energy 
density is obtained as 


1 


(47r + A)<a 


*& +* 53 + *02*03 


Co 2 


t 6n 


n 

e 


(28) 
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Using (28) in we have the pressure 


as 


V = 


(V 0 2 2 + V 0 2 3 + V 02 X 03 )^- 6ra+2 - Yin 

(An + \)V 0 H 2 


l(X 2 2 + X 2 3 + Xo 2 X 03 )t- 6n+2 - V 2 n 


(47r + A )Vg 2 a:i 2 


(29) 


With the help of (27)-(29) from (|T9|, the cosmological constant A is obtained as 

-4 


A 


(47r + A )V£t en + 2 _ 


. X w. . v U ,2,6n + 2 (X 2 2 + X 2 3 +X 02 X 0 s)t-^+ 2 -V 2 n 

(27T + A)(4tt + A)C 0 t - J -( 47r T A Wed? - 


- (tt + ^)(V^ 2 + AT 0 2 3 + X 02 X 03 )t 2 + l(3An + A + 12n7r)C 0 2 nf 6n 


(30) 



t 


Fig. 1. Variation of energy density p against 
time t for A = 1 , a = —0.1 ,Vo = 
1 ,X 02 = 0.01 ,AAi 3 = 0.01 and different 

n(1.2,1.4,1.6,1.8) 



t 


Fig. 2. Variation of Pressure p against 
time t for A = 1 , a = —0.1 ,Vo = 
l,Vo 2 = 0.01,Vo3 — 0.01 and different 
n(1.2,1.4,1.6,1.8) 



Fig. 3. Variation of Cosmological constant A against time £ for A = 1 , a = — 0.1 ,Vo = 1,V 0 2 = 0.01,V O 3 = 0.01 and different 
n(1.2,1.4,1.6,1.8) 


Fig|l]and Fig(2]show the variation of energy density p and pressure p against time t for different values as in the 
figures. Here we noticed that 0 when t oo. In the increase of 77, energy density and pressure increases and 

decreases respectively. Figj3] represents the variation of cosmological constant A against time t for different values as 
in the figures. It is observed that Cosmological constant A is decreasing with the increase of n and with the evolution 
of time it approaches towards zero. 
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4.2 Exponential law Model 


In this case, with the help of (24) in (25), we have found the metric potential as 


Xi(t) = XioV^exp 


-/? Vot + 3Xo*e 
3/3Vb 


_. p -pt 


, i = 1,2,3 


(31) 


The directional Hubble parameters are obtained as H t = | + ■ i = 1,2,3. The Hubble parameter (H), deceleration 

parameter ( q) , expansion scalar (6>) and Shear scalar (a 2 ) are as follows 


rr P , , p o 2 X 2 2 +X 2 3 +X 02 X 03 

H=-,q = -l,0 = P,a = -- 


Here we noticed that a 2 die out for larger values of t. From (31) and (26), the energy density is expressed as 


P 2 = 


{X 2 2 + X 2 3 + X 02 X 03 )e- 2 P t 


V 2 a(4n + A) 


(32) 


Using (32) in & the pressure is expressed 


as 


P = 


(47t + \)Vq 



F 0 2 a(47r + A) 


Vo 2 Vo 3 )e _2/3t 


(33) 


With the help of (31)-(33) from (19), the cosmological constant A is obtained as 




—— X X r\ I ^ 


(34) 



t 


Fig. 4. Variation of energy density p against 
time t for A — 1 , a — 0.1 ,Vo — 

1,X 0 2 — 0.01,V O 3 = 0.01 and different 
£(0.2,0.4,0.6,0.8,1) 



Fig. 5. Variation of Pressure p against time t 
for A = 1 , a = 0.1 ,V 0 = 1,X 02 = 0.01,V 03 = 
0.01 and different £(0.2, 0.4, 0.6, 0.8,1) 
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Fig. 6. Variation of Cosmological constant A against time t for A = 1 , a = 0.1 ,Vo = l,Vo 2 = 0.01,A4)3 = 0.01 and different 
£(0.4,0.6,0.8,1,1.2) 


Fig (4] and Fig (5] show the variation of energy density p and pressure p against time t for different values as in 
the figures. Here we noticed that —>> 0 when t 00 . In the increase of £, energy density and pressure decreases 
respectively. Figj 6 ] represents the variation of cosmological constant A against time t for different values as in the 
figures. It is observed that Cosmological constant A is not approaching towards zero with the evolution of time and 
also it takes negative values. 


5 Field equations and cosmological model for f(R,T ) = fi(R) + f 2 {T) 

In f(R,T) theory, the gravitational field equations ( |15| ) for the choice of fi(R) = A R and / 2 (T) = AT, along with 
Cosmological constant A are given as 


Gij 


8 tt + A 
A 


Rj + 


p-p + 2 A\ 


I 9ij ) 


In this case, the field equations are given by 


Ms V 2 X 3 

VX 2 m 3 + X 2 V 3 


167T + 3AA p 

-^) P+ 2 - A 


X 2 X 3 X 2 X 3 
x 2 + x 3 + X 2 X 3 

X1X3 X1X3 
Xi x 3 + x,x 3 


Xi 

X 1 


V 2 

V 2 


W -2 

X!X 2 


I 67 r + 3A 
2 A 

I 67 r + 3A 
2 A 

I 67 r + 3A 
2 A 


>- 2 - A 


P -A 




(35) 

(36) 

(37) 

(38) 

(39) 


5.1 Power Law Model 


Following the same procedure as in subsection (4.1), we have obtained the same metric potential as in equation (27) 


and the other parameters like energy density p, pressure p and cosmological constant A are expressed as follows 


p 2 = 


2A 


C 0 2 a(87r + A) 


Xq 3 + Af 0 2 2 + Xq 2 Xq 3 Vt 


V 2 n 


t 6 n 


t 2 


(40) 


2A 

rx 0 2 3 + x 2 2 + x 02 Xo 3 

V 2 n] 

J 2A 

[X 2 3 + X 2 2 + X 02 X 03 

V 2 n] 

y 2 (87r + A) 

t 6n 

t 2 _ 

V KM87T + A) 

t 6n 

t 2 _ 


(41) 
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A = 


41^(47r + X)V8^TXV^V 0 t~ 3n+1 - 2i 2 (X 2 3 + x[ 2 + X^X^) 


2£ 2 A(8tt + A)aF 0 2 

+ 4(4t r + A)(X 0 2 3 + Xl 2 + Xo 2 X Q 3 )Xat ~ 6n+2 - 2n(3An + A + 24ri7r)V 0 2 a;A 


(42) 



Fig. 7. Variation of energy density p against 
time t for A = 1 , a m —0.1 ,Vo — 
l,Xo 2 = 0.01,Vo3 = 0.01 and different 
n(1.2,1.4,1.6,1.8) 



t 

Fig. 8. Variation of Pressure p against 
time t for A « 1 , a = —0.1 ,Vo — 
l,Vo 2 = 0.01,Vo3 = 0.01 and different 
n(1.2,1.4,1.6,1.8) 



Fig. 9. Variation of Cosmological constant A against time £ for A = 1 , a = — 0.1 ,Vo = l,-Xo 2 = 0.01,Xo3 = 0.01 and different 
71(1.2,1.4,1.6,1.8) 


Here also we have noticed the similar results as in (4.1). 


5.2 Exponential law Model 


Following the same procedure as in subsection (4.2), we have obtained the same metric potential as in equation (32) 


and the other parameters like energy density p, pressure p and cosmological constant A are as follows 


P = 


1 2(X 2 3 + X 2 2 + Xo 2 Xo 3 )X 


VoePt 


a(87r + A) 


(43) 


P = 


2(V 


03 


x 2 02 


■ X 02 X 03 ) X 


l2(X ( 


Xq2 


' X02X03JX 


V 0 2 e 2 ^( 87 r + A) 


V 0 e& 


a( 87 r + A) 


(44) 
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A — — 


2V2(4tt + A) /2(ATq 3 + AT 2 2 + X 02 X 03 )\e~W 2(X 2 3 + X 2 2 + AT 02 V 03 )(47r + A)e“ 2 ^ 


C 0 2 a(87r + A) 


Vq (87t + A) 


3 


(45) 



t 


Fig. 10. Variation of energy density p 
against time t for A = 1 , a = 0.1 ,Vo = 
l,Xo 2 = 0.01,Vo3 = 0.01 and different 
£(0.2,0.4,0.6,0.8,1,1.2) 


x 10' 3 



Fig. 11. Variation of Pressure p against 
time t for A — 1 , a — 0.1 ,Vo — 

l,Vo 2 = 0.01,Vo3 = 0.01 and different 
£(0.2, 0.4, 0.6, 0.8,1,1.2) 



Fig. 12. Variation of Cosmological constant A against time t for A = 1 , a = 0.1 ,Vo = l,Vo 2 = 0.01,X 03 = 0.01 and different 
£(0.4,0.6,0.8,1,1.2) 


Here also we have noticed the similar results as in (4.2). 


6 Concluding Remarks 

In this article we have the Bianchi type-I cosmological model in f(R,T ) modified gravity for two different classes 
of /(i?, T) in presence of cosmological constant and quadratic equation of state. Here we have discussed two models 
based the expansion law. From both the models, case of /(i?, T) = i?+ 2/(T), we have concluded the following points. 

— In both the models, energy density p is decreasing function of t and p approaches towards zero with the evolution 
of time. 

— In both the models, pressure p is negative and approaches towards zero with the evolution of time. 

— In both the models, cosmological constant A is negative but here we notice that incase of Power law A approaches 
towards zero with the evolution of time where as it does not approaches towards zero with the evolution of time 
incase of exponential law. 
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Similar observation also noticed for the case of f(R, T) = fi(R) + / 2 (T). Here all the observation are in fare agreement 
with the observational data. 
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